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Background One-nucleon removal reactions at or above the Fermi energy are important tools to explore the single-particle 
structure of exotic nuclei. Experimental data must be compared with calculations to extract structure information, 
evaluate correlation effects in nuclei or determine reaction rates for nuclear astrophysics. However, there is insufficient 
knowledge to calculate accurately the cross sections for these reactions. 

Purpose We evaluate the contributions of the final state interaction (FSI) and of the medium modifications of the nucleon- 
nucleon interactions and obtain the shapes and magnitudes of momentum distributions. Such effects have been often 
neglected in the literature. 

Method Calculations for reactions at energies 35 - 1000 MeV/nucleon are reported and compared to published data. For 
consistency, the state-of-the-art eikonal method for stripping and diffraction dissociation is used. 

Results We find that the two effects are important and their relative contributions vary with the energy and with the atomic 
and mass number of the projectile involved. 

Conclusions These two often neglected effects modify considerably the one-nucleon removal cross sections. As expected, the 
effect are largest at lower energies, around 50 MeV/nucleon and on heavy targets. 



PACS numbers: 25.60. Gc,21.65.-f,21. 10. Jx 

I. INTRODUCTION 

Nuclcon knockout reactions in nuclear collisions at and 
above the Fermi energy in nuclei have become an impor- 
tant tool to determine the occupancy of single-particle 
states and the correlation effects in the nuclear many- 
body system (see, e.g. Refs. [H-IH). In peripheral, 
sudden collisions of fast-moving projectiles with a tar- 
get nucleus, a single nuclcon is removed from the projec- 
tile, producing projectile- like residues in the exit channel, 
which are measured. Referred to the center-of-mass sys- 
tem of the projectile, the transferred momentum is k c . 
For the knockout reactions in the sudden approximation, 
this must equal the momentum of the struck nucleon be- 
fore the collision. The standard reaction models assume 
that the ground state of the projectile of spin and parity 
J* can be approximated by a superposition of configu- 
rations of the form [JJ C <E> nlj] J where T£ e denote the 
core states and the nlj are the quantum numbers for the 
single particle wave functions in a spherical mean field 
potential. The measured partial cross sections to indi- 
vidual final levels c of the core allow to extract, by com- 
parison with theoretical calculations, the spectroscopic 
factors for the individual core-single particle configura- 
tions. In complete analogy to the use of angular dis- 
tributions in transfer reactions, the orbital angular mo- 
mentum I is revealed by the shape of the momentum 
distributions -P(k c ). It is obvious then that the accuracy 
of the extracted spectroscopic factors, a measure of the 
occupancy of the single-particle orbitals in nuclei, and 
the conclusions that may follow from these spectroscopic 
factors about correlations inside nuclei depend in direct 
measure on the accuracy of the cross section calculations. 
Similarly, the cross sections for one-proton removal reac- 



tions X — >• Y + p are directly related to the non-resonant 
part of the astrophysical S-factors for the inverse radia- 
tive proton capture Y(p,7)X (sec, e.g., Refs. 
Again, the results and their reliability depend directly on 
the reliability and accuracy of the reaction calculations 
(as discussed in [l(il|). 

The one-nucleon removal cross section is calculated in 
most reaction models as an incoherent sum of the contri- 
butions of all core-single particle configurations making 
the ground state of the fast moving projectile: 

cr_i„ = ^2S(c; nlj)o- sp (nlj), (1) 

where S(c; nlj) and a sp are the spectroscopic factors of 
each configuration and the single particle removal cross 
section, respectively Q. A similar relation is valid for 
momentum distributions. Systematic studies of projec- 
tiles and reactions allow the determination of the order- 
ing, spacing and the occupancy of orbitals, essential in 
assessing how nuclei evolve in the presence of large neu- 
tron or proton excess. Much was done in this respect 
in the last decade in various laboratories. This informa- 
tion can be compared to many-body nuclear structure 
calculations which are now able to reproduce the mea- 
sured masses, charge radii and low- lying excited states of 
a large number of nuclei. It was found that, e.g., for very 
exotic nuclei, the small additional stability that comes 
with the filling of a particular orbital can have profound 
effects upon their existence as bound systems, their life- 
time and structure, and lead to the discovery of magic 
numbers that do not manifest along the valley of stabil- 
ity. 

Extensions of the nucleon knockout formalism includ- 
ing the treatment of final-state interactions have been 
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discussed in Ref. where it is shown that Coulomb 
final-state interactions are of relevance. In the mean- 
time, inclusion of higher-order effects [H, EH and a the- 
ory for two-nucleon knockout flil - Eiij have been devel- 
oped. Knockout reactions represent a particular case for 
which higher projectile energies allow a simpler theoreti- 
cal treatment of the reaction mechanism, due to the sim- 
plicity of the reaction mechanism and the assumption of 
a single-step process. 



A microscopic approach to direct reactions uses an ef- 
fective nucleon-nucleon (NN) interaction (e.g. those of 
Ref. [It]]) to start with. This interaction is often used 
to construct an optical potential with its imaginary part 
assumed to relate to the real part and its strength ad- 
justed to reproduce experimental data. The real and 
imaginary parts of the potential can also be indepen- 
dent as in Refs. 0,0, where the procedure starts from a 
NN effective interaction with independent real and imag- 
inary parts. For collisions at high energies (E > 100), 
it is possible to show that instead of nucleon-nucleon 
interactions one can use nucleon-nucleon cross sections 
as the microscopic input [l8|. In this case, an effective 
treatment of Pauli-blocking of nucleon-nucleon scatter- 
ing is needed, as it manifests through medium modifica- 
tion of nucleon-nucleon cross sections. It is well known 
that medium modification of the nucleon-nucleon cross 
sections is necessary for an adequate numerical model- 
ing of heavy-ion collision dynamics in central collisions. 
In these collisions, the ultimate purpose is to extract in- 
formation about the nuclear equation of state (EOS) by 
studying global collective variables describing the colli- 
sion process. In direct reactions, such as one-nuclcon 
removal reactions, medium effects of NN scattering are 
smaller because mostly low nuclear densities are probed. 
A first study of this effect in knockout reactions was car- 
ried out in Ref. [l9j . Nonetheless, no comparison with 
experimental data was provided. In this work we ex- 
plore further consequences of medium corrections and fi- 
nal state interactions in knockout reactions. We study 
medium effects in the NN cross section in knockout re- 
actions using the methods reported in Ref. (l9| . namely 
with a geometrical treatment of Pauli-blocking and with 
the Dirac-Brueckner theory in terms of baryon densities. 
We also explore the effect of final state interaction, in 
particular the effects of Coulomb distortion in the en- 
trance and final reaction channels. This is of relevance 
as an increasing number of experiments use heavy targets 
with a large nuclear charge. We compare our results of 
knockout cross section and momentum distribution cal- 
culations to a large number of published experimental 
data. The purpose is to improve the accuracy of the ex- 
tracted spectroscopic factors that will lead to better un- 
derstanding nuclear structure and to check and improve 
the reliability of the use of one-nucleon removal reactions 
as indirect methods in nuclear astrophysics. 



II. MEDIUM AND DISTORTION EFFECTS 

The geometrical treatment of Pauli corrections is per- 
formed using the isotropic NN scattering approximation 
because the numerical calculations can be largely sim- 
plified if we assume that the free nucleon-nucleon cross 
section is isotropic. In this case, a formula which fits the 
numerical integration of the geometrical model reads (l9j 

a NN (E,p p ,p t ) = <J f ^ e (E) 
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where E is the laboratory energy in MeV, p — (p p + 
Pt)/po- P< = min(p p , p t ), p i=p .t is the local density of 
nucleus i, and po = 0.17 fm~ 3 . The parameters and 
models for the p p and p t densities which have been used 
to describe the nuclei in this work are presented in Ta- 
ble Ell 

The Brucckner method goes beyond the simple ge- 
ometrical treatment of Pauli blocking. Some of the 
Brueckner results that we used in this analysis have been 
reported in Refs. 0, Hl|, where a simple parametriza- 
tion was given. It reads (the misprinted factor 0.0256 in 
Ref. [U has been corrected to 0.00256) 
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The limits of validity of this parametrization are clearly 
associated with the limits of validity of the Brucckner cal- 
culations, which are valid only below the pion-production 
threshold. A modification of this parametrization was in- 
troduced in Ref. [22[ and consists in combining the free 
nucleon-nucleon cross sections parametrized in Ref. (23j 
with the results of Brueckner theory reported in Refs. 

mm. 

Current theoretical models for the calculation of mo- 
mentum distributions and cross sections in high-energy 
nucleon-removal reactions follow a semiclassic al p roba- 
bilistic approach, described, e.g., in Refs. 0, [Hf . The 
method relies on the use of "survival amplitudes" (or S- 
matrices) in the eikonal approximation, 



Si(b) = exp[«x(6)] = cxp 



i 

hv 



U iT (r)dz 



(4) 



where r = \Jb 2 + z 2 , and UjT is the particle(i)-target(T) 
optical potential. In Ref. 18[, a relation has been de- 
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FIG. 1: (Color online). Top and right scale: Ratio between 
average in-medium and the free nucleon-nucleon cross section 
as a function of the impact parameter. Dashed and dotted 
curves are for core-target and valence nucleon-target average 
nucleon-nucleon cross sections, respectively. Notice that the 
target center of mass is located on the right of the top axis 
scale. Bottom and left scale: Radial wave functions in arbi- 
trary units (solid curves) for the valence nucleon-core system 
and for a few representative reactions considered in this work. 
We have taken only one configuration in cases of systems with 
multiple configurations. 



veloped between the optical potential and the nucleon- 
nucleon scattering amplitude. Such a relation is often 
referred in the literature as the "t-pp approximation" . 
The t-pp approximation is the basis of most calculations 
of elastic and inelastic scattering involving radioactive 
nuclei, as experimentally deduced optical potentials are 
not often available. In this approximation, the eikonal 
phase becomes 



X(b) 



1 



dq q p P (q) Pt (q) Inn (q) Jo (qb) , (5) 



where p Pjt (q) is the Fourier transform of the nuclear den- 
sities of the projectile and target, and /at n (<?) is the high- 
energy nucleon-nucleon scattering amplitude at forward 
angles, which can be parametrized as 



/nn (q) 
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cussed in Refs. [ll|,|26[. Besides Coulomb repulsion, in- 
cluded by modifying the straight-line trajectories accord- 
ingly, we have also modified the integral in Eq. (j4]) by us- 
ing the optical potential including the Coulomb potential 
that modifies the S-matrices according to Si = S^ ■ Sf . 
The Coulomb phase in Sf is calculated by assuming a 
uniform charge distribution with radius R, and is given 
by 

X c(b) = 2 V {e(b-R) ln(kb) 

+ Q(R - b) \\n(kR) + ln(l + y/l - b 2 /R 2 ) 

-^l-V/R 2 -^{l-b 2 /R 2 f/ 2 | , (7) 

where is the step function. The value R is chosen to be 
small enough so that the nuclear S-matrices are basically 
zero below b = R because of the strong absorption at 
small impact parameters. 



III. RESULTS AND DISCUSSION 

In this section the results for momentum distributions 
and nucleon-removal cross sections arc compared to sev- 
eral experimental data. The focus is the importance of 
medium corrections of nucleon-nucleon cross sections and 
of Coulomb distortions. Both effects are expected to de- 
crease as the bombarding energy increases. It is impor- 
tant to include such effects in order to minimize the un- 
certainty in the extraction of spectroscopic factors, espe- 
cially at low bombarding energies. To substantiate this 
assertion, we analyze low energy data on knockout re- 
actions and compare to high energy data. In order to 
identify the separate contribution of these two factors, 
we do not vary the geometry of the nucleon binding po- 
tentials used to calculate the single particle radial wave 
functions. That was identified in the literature as an- 
other major factor in the calculations, leading to large 
variations in the extracted spectroscopic factors. 

The relevance of medium corrections arc moti- 
vated by the effect summarized in Figure [TJ The 
cases are shown in Figure [1] are for 12 C( 17 C, 16 B) at 
35 MeV/u, 9 Be( 15 0, 14 N) at 56 MeV/u, 12 C( 23 Al, 22 Mg) 
at 50 MeV/u and 9 Be( n Be, 10 Be) at 60 MeV/u. Later we 
will discuss more results for each system separately. The 
dashed and dotted curves show the ratio between the av- 
erage in-medium and the free nucleon-nucleon cross sec- 
tion as a function of impact parameter (see top scale). 
We define the average nucleon-nucleon cross section at 
the distance of closest approach between the projectile 
and the target using Eq. and the definition 



ctnn (i + a NN )cxp(-f3 NN q 2 ) . (6) 



a(E,b)) 



J d 3 r p p (v)pt(v - b) o- NN (E, p p , p t ) 
J d 3 rp p (r)p t (r - b) 



(8) 



There are many ways of introducing final state inter- 
actions in direct nuclear reactions, some of which are dis- 



wherc b is the impact parameter vector, perpendicular 
to the beam axis. 
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Further, as has been shown in Ref. the transverse 

momentum distributions are strongly affected by both 
nuclear and Coulomb elastic scattering. For heavier tar- 
gets the distortions are predominantly due to Coulomb 
repulsion 
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FIG. 2: Upper panel: Display of Coulomb scattering ef- 
fects in longitudinal momentum distributions for the reaction 
T( 17 C, 16 B)X at 35 MeV/u as a function of target T. The solid 
(dashed) curves are calculations with (without) Coulomb dis- 
tortion. Lower panel: Same as above, but for 12 C target and 
for different beam energies. 



In Figure [TJ the dashed and dotted curves are for 
core-target and valence nucleon-target average nucleon- 
nucleon cross sections, respectively. Notice that the tar- 
get center of mass is located on the right of the top axis 
scale. Also shown in the figure are the radial wave func- 
tions (solid curves) for the valence nucleon-core system 
and for a few representative reactions considered in this 
work. For simplicity, in this figure we have used only 
one of the main configurations for the projectile ground 
state (more detailed and complete calculations are re- 
ported later on in this section). The binding energies 
"effectively" decrease from the top to low panels. We 
mention "effectively" because, although the binding en- 
ergy of the valence proton in 23 Al is smaller than for the 
valence neutron in 11 Be, the Coulomb barrier creates an 
effectively larger binding in 23 Al. 

It is clearly noticeable from FigureQ]that the wavefunc- 
tions of weakly bound systems extend far within the tar- 
get where the nucleon-nucleon cross sections are strongly 
modified by the medium. Momentum distributions and 
nucleon removal cross sections in knockout reactions are 
thus expected to change appreciably with the inclusion 
of medium corrections of nucleon-nucleon cross section. 
Such corrections are also expected to play a more signif- 
icant role for loosely-bound systems. 

In the following, we discuss Coulomb corrections. Here 
we consider the simplest and most straightforward cor- 
rection one can do, namely the inclusion of a Coulomb 
phase, which accounts for the distortion of the clastic 
scattering of the core fragment. It has been usually taken 
for granted that longitudinal momentum distributions 
are little affected by elastic scattering of the core frag- 
ment, the reason being that the longitudinal forces acting 
on the core fragment reverse sign as the projectile passes 
by the target, leading to a reduced distortion effect [l[. 



It is worthwhile mentioning that the im- 
plications of the findings on Coulomb distortion effects 
presented in Ref. [ll| have been neglected in the liter- 
ature. In order to avoid dealing with the effects of the 
Coulomb scattering, experiments are usually performed 
with light targets, such as 9 Be and relatively high ener- 
gies, E > 50 MeV/nucleon. In this work we show that 
these argumentations are not always valid and need to 
be studied with care. 

As discussed in the previous section, in the presence of 
the Coulomb field the eikonal S-matrices factorize as the 
product of the nuclear and the Coulomb contributions: 
S(b) = S n (b)Sc(b). Although this does not make any dif- 
ference for the total stripping cross sections (see Eq. (20) 
of Ref. [lj| ) , it has an impact on the diffraction dissocia- 
tion cross section (through the second term of Eq. (21) of 
Ref. [lj|). This means that not only transverse, but also 
longitudinal momentum distributions will be affected by 
the Coulomb field. This is shown in Figure [2] for the 
longitudinal momentum distributions of several systems 
which we will consider in details later in this section. It 
is evident from the upper panels of this figure that lon- 
gitudinal momentum distributions in knockout reactions 
T( 17 C, 16 B)X (and their total cross sections) are strongly 
influenced by the Coulomb field of the target T at bom- 
barding energies of 35 MeV/nucleon. The solid (dashed) 
curves are calculations with (without) the inclusion of 
Coulomb scattering. It is also evident that even for the 
case of light targets, such as 9 Be and 7 Li, the distribu- 
tions change appreciably. The lower panels show calcu- 
lations for the same reaction, but for 12 C targets and as 
a function of the bombarding energy It is clear that dis- 
tortions are important even for usually considered "safe" 
energies, such as 100 MeV/nucleon. 

We found that the effect of Coulomb scattering is rela- 
tively larger for systems with smaller sizes. This is illus- 
trated in Figure [3] where we present our calculations for 
the total nucleon removal cross sections for the reactions 
12 C( 17 C, 16 B) (solid curve) and 12 C( 23 Al, 22 Mg) (dashed 
curve) at 35 MeV/u. We artificially vary the binding- 
energy of the proton in 17 C and in 23 Al. As the sepa- 
ration energy increases so does the percent difference of 
the cross section, (a — cr no coui)/c, where a no coui is a 
the cross section without the Coulomb scattering phases. 
With increasing separation energy the relative valence 
nucleon-core distance decreases and the nucleon removal 
cross section decreases, but the relative importance of 
the Coulomb scattering increases. For very small ener- 
gies the effects of Coulomb dissociation (not considered 
here) should also become relevant and increase the mag- 
nitude of the cross sections. The relative importance of 
the Coulomb scattering for the removal cross sections de- 
creases with the bombarding energy. This is shown in 
Figure [3] with the calculation for 12 C( 23 Al, 22 Mg) at 50 



FIG. 3: Total nucleon removal cross sections for the reac- 
tions 12 C( 17 C, 16 B) (solid curve) and 12 C( 23 Al, 22 Mg) (dashed 
curve) at 35 MeV/u. We artificially vary the separation en- 
ergy S of the proton in 17 C and in 23 Al. The dotted curve 
shows the calculation for 12 C( 23 Al, 22 Mg) at 50 MeV/u. 



McV/u (dotted curve). 

These preliminary discussions support our conclusion 
that both medium effects and Coulomb distortion play 
a relevant role in knockout reactions. Next we consider 
a series of published data for which neither medium or 
Coulomb corrections were included. We thus quantify the 
changes in the extracted values of spectroscopic factors in 
case these effects were to be included in the experimental 
analysis. 



A. 12 C( 23 Al, 22 Mg)X at 50 MeV/u 

Recently, the 12 C( 23 Al, 22 Mg)X knockout reaction has 
been studied at 50 MeV/nucleon to investigate the 
ground state properties of 23 Al It was shown that the 
ground-state structure of 23 Al is a configuration mixing 
of a d-orbital valence proton coupled to four core states 
of 22 Mg - 0+ , 2+, 4+, 4+. The ground state spin and 
parity of 23 Al as J 71 ' = 5/2 + has been confirmed. The 
results of Ref. [1] are shown in Fig. [4] and Table Q] This 
experiment had the advantage that exclusive measure- 
ments were done and momentum distributions were de- 
termined for the four major configurations in the ground 
state of the projectile ( 23 A1). 

In this work, we have analyzed the 12 C( 23 Al, 22 Mg)X 
system to check the relevance of Coulomb and medium 
effects. The wave functions for the valence proton 

were generated in a spherical Woods-Saxon (WS) poten- 
tial with the parameters given in Table |TTJ 

In the optical limit of the Glauber theory and the t- 
pp approximation (explained in detail in Refs. (l7l. Il8j]). 
the eikonal phase is obtained from the input of the nu- 
clear ground state densities and the energy dependent 
nucleon-nucleon cross sections. The ground state den- 
sity parameters and models used in this work are shown 




p„ (MeV/c) 

FIG. 4: (Color online). Exclusive longitudinal momentum 
distributions. Shaded areas correspond to lcr deviation in the 
spectroscopic amplitudes 0]. 



in Table [TTTl and our results are presented in Fig. [5] and 
Table H 

To understand the effects of medium and Coulomb cor- 
rections, we have performed the calculations with differ- 
ent inputs. We show in Figure [5] the calculations with 
both Coulomb and medium corrections (solid curve), cal- 
culations without any medium corrections (dashed lines) , 
calculations that exclude Coulomb distortions but keep 
medium corrections (dashed-dotted curve), and calcu- 
lations without either Coulomb or medium corrections 
(dotted curve). 

The numerical results for the single particle cross sec- 
tions with different configurations are shown in Table HI 
For each of the four configurations - 22 Mg - 0+ s , 2+, 4+, 
4% - the corresponding relative differences between full 
calculations and calculations without Coulomb correc- 
tions are found to be 14%, 16%, 18% and 19%, respec- 
tively, whereas between full calculations and calculations 
without medium corrections the corresponding percent- 
age differences are found to be 25%, 21%, 17% and 14%, 
respectively. 



B. 9 Be( 15 0, 14 N)X at 56 MeV/u 

One-proton removal reaction from 15 O on a Be tar- 
get has been measured at 56 MeV/nucleon and the total 
knockout cross section is reported as 80±20 mb in Ref. 
[3lj . The authors were able to explain the orbital occu- 
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Configuration 



E x 
fkeVl 



^Mg(0+)®7Tl d 

22 Mg(2+) ® ir ld 
22 Mg(4+)®^ ld5/2 
22 Mg(4+)®^ ld5/2 
inclusive 



[mb] 



tit 
[mb] 



5', 



5/2 
5/2 



18.5 ± 1.2 
1247 39.3 ± 1.2 
3308 9.5 ± 0.9 
5293 10.9 ± 0.9 

78.3 ± 0.4 77.8 ± 0.7 



19.1 ± 2.0 0.36 0.45 ± 0.07 

39.0 ± 2.2 0.92 1.15 ± 0.18 

9.3 ± 0.9 0.27 0.34 ± 0.06 

10.4 ± 1.6 1.43 0.50 ± 0.09 











a sp (nlj) 


[mb] 




Configuration 


[keV] 


Full 


no medium no Coulomb 


Free 


"Mg(0+ ) C 
22 Mg(2+) « 
22 Mg(4+) 6 
22 Mg(4+) 6 


3 7Ti d5/2 
3 7Tld 5/2 
3 Kld 5/2 
3 7T ld5/2 




1247 
3308 
5293 


27.1 
23.7 
20.4 
18.4 


33.8 
28.7 
23.9 
21.0 


23.2 
19.9 
16.7 
14.8 


30.1 
25.1 
20.4 
17.6 



TABLE I: Upper table 
separately. 



Previous work [Jj]. Lower table - Present work. Single particle cross sections are shown for each case 



J 7T 


V 


Ro 


a 


v s0 


Rao 




R c 






(MeV) 


(fin) 


(fin) 


(MeV) 


(fin) 


(fin) 


(fin) 


(MeV) 


\ w Be{r)®v2s 1/2 ) 




61.13 


2.70 


0.52 








2.70 


0.504 


\ 1/l N{J 7T )(gmnlj) 




48.36 


2.94 


0.60 








2.94 


7.297 


(1P3/2) 


48.36 


2.94 


0.60 








2.94 


7.297 


| 1D B(J")®^lp3/ 2 } 




79.46 


2.80 


0.50 


35.0 


2.80 


0.50 


2.80 


23.330 




80.35 


2.80 


0.50 


35.0 


2.80 


0.50 


2.80 


23.979 




80.75 


2.80 


0.50 


35.0 


2.80 


0.50 


2.80 


24.273 




81.85 


2.80 


0.50 


35.0 


2.80 


0.50 


2.80 


25.078 




82.17 


2.80 


0.50 


35.0 


2.80 


0.50 


2.80 


25.318 


3 2 " 


79.93 


2.85 


0.50 


25.0 


2.80 


0.50 


2.80 


26.066 


|^Mg(J' r )®7rld 5/2 ) 




54.60 


3.36 


0.60 


5.0 


3.36 


0.60 


3.36 


0.141 


2+ 


56.96 


3.36 


0.60 


5.0 


3.36 


0.60 


3.36 


1.388 


jr 


60.67 


3.36 


0.60 


5.0 


3.36 


0.60 


3.36 


3.449 


4 2 + 


64.07 


3.36 


0.60 


5.0 


3.36 


0.60 


3.36 


5.434 




1/2/ , 
' (9.3) 


42.40 


3.65 


0.70 








3.65 


3.610 


| 3a Mg( J^i/nZj) 


or 9 . s) ( irf 3/ 2 ) 
















2.21 


3 _ (2p 3 /a) 


79.92 


3.35 


0.70 


10.0 


3.30 


0.70 


3.35 


5.07 


3"(l/7/a) 


86.63 


3.35 


0.70 


10.0 


3.30 


0.70 


3.35 


5.07 


2j(2s 1/2 ) 


51.55 


3.35 


0.70 


10.0 


3.30 


0.70 


3.35 


5.22 



TABLE II: Bound state potential parameters for the systems studied in the present work. S e ff is effective separation energy: 
S e ff=Si+Ex° re where i — proton or neutron and E% or£ is core excitation energy. 



pancy of valence protons with a pure lpi i-x single particle 
state using a Glauber reaction model. Their calculations 
imply that the 1^3/2 state could also have a small contri- 
bution because the calculations with only the lpi / 2 state 
yield a narrower momentum distribution than observed 
in the experiment. The physical implication of this is a 
possible knockout from more deeply bound protons in the 
I/J3/2 state. The contributions from each of the p-states 
yield spectroscopic factors of 1.27(9) and 0.100(75) for 



the lpi/2 and the IP3/2 orbital, respectively (Ref. [3l| 
and references therein) . 

We have followed the interpretation of Ref. [3l| and 
calculated the one-proton removal cross sections for the 
same reaction with the same orbital occupancy assump- 
tion. The parameters are shown in Tables [TT| and IIIII 
Our calculations with both Coulomb and medium cor- 
rections by slightly changing the spectroscopic factors as 
1.42 and 0.13 are in agreement with the results of Ref. 
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-200 200 
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FIG. 5: (Color online). Comparison of experimental data of 
figure [4] and calculations for exclusive longitudinal momen- 
tum distributions in the knockout reaction 12 C( 23 Al, 22 Mg)X 
at 50 MeV/nucleon. The solid line has both Coulomb and 
medium corrections. The dashed-curve has no medium cor- 
rections. The dashed-dotted line includes calculations with- 
out Coulomb corrections. The dotted curve neither includes 
medium effects nor Coulomb corrections. 







(fin) 


(fm) 


(fm) 


(Target) 










9 Be 


HO 


2.50(9) 


1.77(6) 


0.631 


12 C 


MHO 


2.440(25) 


1.635(9) 


1.403(16) 


(Projectile) 










14 N 


HO 


2.540(20) 


1.729(6) 


1.291 


16 B 


LDM 








22 Mg 


HFB [27-29J 


2.92 






23 o 


LDM 








32 Mg 


2pF 


3.06(5) 


3.05(5) 


0.523(12) 



TABLE III: Ground state densities adapted from elec- 
tron scattering data in Ref. [3(J| where c and z are for 
the 2pF model, a and a for the HO and MHO models. 
22 Mg ground state density obtained with a Hartree-Fock- 
Bogoliubov (HFB) calculation [27|-[2?|. - LDM: Liquid Drop 
Model. - HO: Harmonic-oscillator model. - MHO: Modified 
harmonic-oscillator model. - 2pF: Two-parameter Fermi dis- 
tribution. 
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[3l| . The calculated one-proton removal cross sections 
are 74.10 mb, 70.26 mb, 87.68 mb and 84.26 mb with both 
Coulomb and medium corrections, no Coulomb correc- 
tions, no medium corrections and neither medium effects 
nor Coulomb corrections, respectively. The difference be- 
tween full calculations, including medium and Coulomb 
scattering effects, and calculations without Coulomb cor- 
rections is of the order of 5%, and between full calcula- 
tions and calculations without medium effects is nearly 
18%. This is remarkable even though it fits again within 
the error of the total knockout cross section experimen- 
tal data. We thus conclude that for this case, medium 
effects and Coulomb distortion do not have a sizable im- 
pact on the extraction of spectroscopic factors. However, 
one can easily see from Fig. [5] that the data shows an 
asymmetry which can only be explained with inclusion 
of higher-order effects in the reaction mechanism. Dis- 
tortions will be manifest due to continuum-continuum 
coupling of states involving the interaction of core with 
the valence proton. These mechanisms have not been 
considered in the present work. 



FIG. 6: (Color online) Longitudinal momentum distributions 
for the 9 Be( ls O, 14 N)X reaction at 56 MeV/u. Solid lines rep- 
resent calculations include both Coulomb and medium correc- 
tions. Dashed lines stem from calculations that do not include 
medium corrections. Calculations denoted by dashed-dotted 
curves are performed without Coulomb corrections. The dot- 
ted curve does not include medium effects or Coulomb cor- 
rections. The data are taken from Ref. [3lj |. 



C. 12 C( 17 C, 16 B)X at 35 MeV/u 

1. Transverse momentum distributions 

One-proton removal reaction from 17 C, with a separa- 
tion energy of 23 MeV, has been measured in the reac- 
tion 12 C( 17 C, 16 B)X at 35 MeV/nucleon with the goal to 
understand the low-lying structure of the unbound 16 B 
nucleus. Using this reaction, Ref. [32j studied the un- 
bound B+n system with the assumption of a e?-wavc 
neutron decay. Our interest is to compute the transverse 
momentum distribution of the 16 B fragment following the 
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and without Coulomb distortion. The reason for this 
difference is that the Coulomb distortion and repulsion 
effectively increases the collision distance at the small 
impact parameters needed to remove a strongly bound 
nucleon. This was not observed in the previous case 
( 9 Be( 15 0, 14 N)X at 56 MeV/u) because of the small nu- 
clear binding in that case. We have also observed that 
this effect sharply reduces the calculated cross sections 
and the removal is more effective as the bombarding en- 
ergy decreases. 



FIG. 7: (Color online). Transverse momentum distributions 
for the 12 C( 17 C, 16 B)X system at 35 MeV/u. Solid lines rep- 
resent calculations including both Coulomb and medium cor- 
rections. Dashed lines stem from calculations that do not 
include medium corrections. Calculations represented by 
dashed-dotted curves are performed without Coulomb cor- 
rections. The dotted curve does not include medium effects 
nor Coulomb corrections. The data are taken from Ref. 1321. 



same assumptions as in Ref. [32| in order to study the 
consequences of medium and Coulomb corrections. The 
configuration of the proton removed from 17 C is assumed 
to be 

| 17 C)=a 1 | 16 B(0-)® 7 rlp 3 / 2 ) 
+ a 2 | 16 B(3r) <8 7rlp 3/2 ) + a 3 | 16 B(2r) ® 7rlp 3/2 ) 



+ a 4 PB(2 2 -) ® nlp 3/2 ) + a 5 \ Lb B{l^) <g> Trlp 3/2 ) 



a 6 | ib B(3 2 -)®^lp 3 /2> 



(9) 



where is the spectroscopic amplitude for a core-single 
particle configuration i = [c® nlj). 

Using spectroscopic factors obtained by means of a 
shell- model calculation with the WBP interaction (33|, 
Ref. [32| obtained a good agreement between data and 
calculated transverse momentum distributions. But the 
measured total cross section is 6.5(1.5) mb against a the- 
oretical result of 24.7 mb. The explanation of this large 
difference is proposed in Ref. [13] as a reduction of the 
spectroscopic factor by 70% for strongly bound nucleon 
systems. After this spectroscopic reduction is accounted 
for, the theoretical estimates for the cross section be- 
comes 7.5 mb, in reasonable accordance with the data. 

In the present work, we do not elaborate on the as- 
sumption introduced in Ref. [32j], and we use the same 
configuration and spectroscopic factors as in (32|. The 
proton binding potential parameters are shown in Ta- 
ble |TT] which are adjusted to obtain the effective separa- 
tion energies. The ground state densities are also listed 
in Table IIIII Here, as it is shown in Figure we find 
that medium corrections change the cross sections by 
5%, but the Coulomb corrections have a very large effect 
which is almost 60% between calculations with Coulomb 



2. Longitudinal momentum distributions 



We have made a more systematic analysis to under- 
stand the reason of the effect discussed in the previous 
subsection. We have observed that the strong depen- 
dence on Coulomb distortions are also present in longitu- 
dinal momentum distributions. It has long been thought 
that longitudinal momentum distributions are free from 
uncertainties related to the knowledge of the optical 
nucleus-nucleus potentials when compared to the trans- 
verse distributions. This was first shown in Ref. Here 
we report calculations for the same | 16 B(0~) (g> irlp 3 / 2 ) 
configuration, with the same parameters and ground 
state densities, as discussed in the previous subsection. 
We find that the longitudinal momentum distributions 
show a similar effect as the transverse momentum distri- 
butions, as can be seen in Fig. [2] 

The source of this difference stems from the diffraction 
dissociation contribution to the cross sections. To sub- 
stantiate our claim, we have looked at the details of the 
knockout cross section which has two parts for the pro- 
duction of a given final state of the residue. The most 
important of the two, commonly referred to as stripping 
or inelastic breakup, represents all events in which the 
removed nucleon reacts with and excites the target from 
its ground state. The second component, called diffrac- 
tive or clastic breakup [35| . represents the dissociation 
of the nucleon from the residue through their two-body 
interactions with the target, each being elastically scat- 
tered. We notice that the total stripping cross section is 
given by [ll| : 



2n 



0~str 



S(c;nlj)——— / db n b n l-\S n {b n )\ 
11 + 1 m Jo 1 



x / d 3 r |S c (& c )| 2 |^ m (r)| 2 , (10) 



whereas the integrated diffraction dissociation cross sec- 
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FIG. 8: (Color online). Longitudinal momentum distribu- 
tions of for the reaction 9 Be( n Be, 10 Be) at 60 MeV/nucleon. 
Solid lines represent calculations include both Coulomb and 
medium corrections. Dashed lines stem from calculations that 
do not include medium corrections. Calculations denoted by 
dashed-dotted curves are performed without Coulomb correc- 
tions. The dotted curve does not include medium effects or 
Coulomb corrections. The data is taken from Ref. 13611. 



tion is given by [26J: 



CTdif = 5(c; nl j) 2l + / db nK 

m 



d 3 r 



S n (b n ) S c (b c )ipi 



E 



d 3 ripi m > (r) S c (b c ) S n (b n )tpi ri 



One can see from these expressions that the stripping 
cross sections are not affected by the Coulomb distortions 
because this distortion is manifest through a real phase in 
the cikonal S-matrices calculated in the Glauber approxi- 
mation. The magnitude of the cross sections are therefore 
not changed, as the square of the S-matrices entering Eq. 
(|10|) are only changed by the imaginary part of the po- 
tential entering Eq. (U). On the other hand, the second 
term of the diffraction dissociation cross sections in Eq. 
(fTTjl is appreciably modified by the Coulomb phase fac- 
tor. As seen from Figure [5J the effect gets smaller with 
decreasing target atomic number because the Coulomb 
phase increases, or when the beam energy increases be- 
cause then the Coulomb recoil becomes irrelevant. 



D. 9 Be( 11 Be, 10 Be)X at 60 MeV/u 

In order to further understand the dependence of the 
Coulomb distortion on the nuclear binding, we consider 
the reaction 9 Be( n Be, 10 Be)X at 60 MeV/u which can 
be modeled by a core plus valence system with the as- 
sumption lo Bc gs (0 + ) + n in 2s\/2 orbital for the ground 



state of n Bc gs (l/2+) (5 n = 0.504 MeV). Here we use the 
same Woods-Saxon potential parameters for the bound 
state as published in Ref. HH]: (i? = 2.70 fm, a = 0.52 
fm) . In Figure [5] and Table IIVI we present our results for 
the the neutron removal longitudinal momentum distri- 
bution of 60 MeV/nucleon n Be projectiles incident on 
9 Be targets. 



= &dif + 0"sfi- 


"C("C/°B) 
Full no Coul no med 


u Be( ii Be, lu Be) 
Full no Coul no med 


Strip. [mb] 
Diff. [mb] 
Total [mb] 


7.56 7.56 5.63 
18.42 2.40 19.15 
25.98 9.96 24.78 


122.5 122.5 164.1 
49.6 43.9 97.3 
172.1 166.4 261.4 



TABLE IV: The cross sections calculated for the systems, 
12 C( 17 C, 16 B) at 35 MeV/nucleon and 9 Be( 11 Be, 10 Be) at 
60 MeV/nucleon. 



It is evident from Fig. [T] that 17 C has the smallest 
"effective" size and that 11 Be has the biggest size among 
the low energy systems in this study. The nuclear size 
is important for low energy cases because the diffraction 
dissociation becomes dominant when the nuclear size is 
smaller, but the stripping becomes dominant when the 
nuclear size is bigger. The reason for this is that a large 
projectile feels the nuclear interaction already at large 
impact parameters. A small projectile can come closer 
to the target where the Coulomb interaction is larger. 
The evidence of this can be seen in Table UVl It is thus 
clear why medium and Coulomb corrections are more 
important in the 9 Bc( 11 Bc, 10 Bc) and the 12 C( 17 C, 16 B) 
cases, respectively. 



E. 12 C( 24 0, 23 0)X at 920 MeV/u 

The momentum distribution of the one-neutron re- 
moval residues from the 12 C( 24 0, 23 0)X reaction was 
measured for the first time at 920 MeV/nucleon and 
reported in Ref. [37j ■ The data could bc explained 
with a spectroscopic factor 5=1. 74(19) of an almost 
pure 2si/2 single-particle state for the valence neutron. 
This work, together with recent theoretical calculations, 
suggests that 24 O is a newly discovered doubly magic 
nucleus. The one-neutron removal cross section was 
found to be 63(7) mb. The calculations in Ref. [37} 
were based on a few-body Glauber formalism [38j for 



(a) 23 
9S (5/2 + ) 



B -(l/2 
- n in Id 



n in 2s x /2 or- 
orbital. The 



two configurations 

bital and (b) O gs \o/ ^ • ) -t- n m ia 5 / 2 
wave functions for the configurations are obtained with 
a Woods-Saxon potential by adjusting the depth of the 
potential to reproduce the one-neutron separation energy 
5„=3.61(27) MeV [13] ■ Using a pure 2si/ 2 configuration 
with 5=1 leads to a cross section of 34 mb. The calcula- 
tion is in agreement with the data when it is multiplied 
by 5=1.74(19). This large spectroscopic factor indicates 
that the single-particle strength of the valence neutron is 
strongly weighted in the 2si/2 state. 
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FIG. 9: (Color online). The longitudinal momentum distribu- 
tions for 12 C( 24 0, 23 0)X reaction at 920 MeV/nucleon. The 
curves are calculated with the free NN cross sections (solid), 
with a geometrical account of Pauli blocking (dashed) , a phe- 
nomenological fit from Ref. [z| (dotted), and a correction 
from Brueckner theory (dashed-dotted) . The data has been 
taken from Ref. 13711 - 



FIG. 10: (Color online). The inclusive longitudinal mo- 
mentum distributions for the 12 C( 33 Mg, 32 Mg)X system at 
898 MeV/nucleon. The data has been taken from R. Kanungo 
et al. [431 . The curves are calculated with the free NN cross 
sections (solid) , with a geometrical account of Pauli blocking 
(dashed), a phenomenological fit from Ref. (dotted), and 
a correction from Brueckner theory (dashed-dotted). 



In the present work we have reproduced the data of 
Ref. [13| also by assuming a 2s\/ 2 orbital only. The po- 
tential parameters for the bound state wave function are 
given in Table [TT] and the ground state density for the 
23 O gs core is obtained using liquid droplet model (LDM) 
densities [39| , as indicated in Table IIIII To understand 
the differences between medium effect models, four dif- 
ferent calculations including Coulomb corrections have 
been made for this system. The calculated one-neutron 
removal cross sections are 59.04 mb, 54.48 mb, 79.06 mb 
and 53.80 mb using free [lj|, Pauli corrected (Eq. 2), 
Brueckner (Eq. 3) and phenomenological parameteriza- 
tions [22j of the nucleon-nucleon cross sections, respec- 
tively. Except for the result obtained with the Brueckner 
theory, they are all in agreement with the previous work 
and with the data. The relative difference between the re- 
sults obtained using Brueckner corrections and with free 
nucleon-nucleon cross sections is about 33%. However, 
we do not consider a real discrepancy, as the Brueck- 
ner parametrization have been extrapolated well beyond 
their validity. Brueckner calculations are limited by the 
pion-production threshold, and should only be valid for 
projectile energies below 300 MeV/nucleon. 

Thus we verify that the experimental data for the re- 
action 12 C( 24 0, 23 0)X at 920 MeV/u is well reproduced 
with the use of free nucleon-nucleon cross sections. The 
changes introduced by Pauli-blocking with the geometric 
model are small, and the phenomenological account of 
medium effects at this high energy also basically agree 
with the results using free cross sections. 



F. 12 C( 33 Mg, 32 Mg)X at 898 MeV/u 

The ground state structure of 33 Mg, a nucleus belong- 
ing to the N = 20 island of inversion, has been studied 
in Ref. [40T | by means of nucleon-removal reactions on 
a carbon target at 898 MeV/nucleon. The longitudinal 
momentum distribution of the 32 Mg core was measured 
and the one-neutron removal cross section was found to 
be 74(4) mb. Most of the contribution to the ground 
state structure of 33 Mg was shown to arise from the 2p 3 / 2 
orbital. 

The longitudinal momentum distribution obtained in 
Ref. cannot be reproduced with a pure single particle 
state. It has been discussed in details in Ref. j4Q| the rea- 
son why a configuration mixing of different single-particle 
states is needed. Two different configuration mixings for 
the ground state of 33 Mg were assumed. The first one is 

| 33 Mg gs (3/2-)) = 

ai| 32 Mg(2+) ® v2p 3/2 ) + a 2 | 32 Mg(l-) ® u2s 1/2 ) 
+ a 3 | 32 Mg(2+) <g> ulf r/2 ) + a 4 | 32 Mg(l-) <g> uld 3/2 ) 
+ a 5 \ 32 Mg(gs) ® v2p 3/2 ) (12) 

and the second is 

| 33 Mg gs (3/2+)) = ai | 32 Mg(3+) ® v2p 3/2 ) 
+ a 2 | 32 Mg(2+) ® v2s y2 ) + a 3 | 32 Mg(3+) ® ulf 7/2 ) 
+ a 4 \ 32 Mg(gs)®vld 3/2 ), (13) 

where Oj are the spectroscopic amplitudes for each single- 
particle orbital. The values of the corresponding spec- 
troscopic factor Si were found by \ 2 minimization and 
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FIG. 11: (Color online). The total reaction cross section 
of the p + 12 C taken from Ref. [4lJ. The curves are cal- 
culated with the free NN cross sections from Ref. [i^l (solid), 
with a geometrical account of Pauli blocking (dashed) , a phe- 
nomenological fit from Ref. [z| (dotted), and a correction 
from Brueckner theory (dashed-dotted). The triangle-dotted 
curve is calculated with the same free NN cross sections from 
Ref. [42[, but with an another HFB calculation [13] for the 
12 C ground state density. 



their values for the second configuration are Si = 2.2 
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In our calculations we have chosen the second config- 
uration set used in Ref. [loj], Eq. [T3l since the 33 Mg 
ground state is usually accepted to be ,V = 3/2 + . We 
apply the same procedure as described before to obtain 
bound state wave functions and cikonal phases. The pa- 
rameters for the bound state potentials and ground state 
densities are shown in Table Ql] and Table IIII) respec- 
tively. We have used nearly the same spectroscopic fac- 
tors within the error bar range of Ref. [4(| to make a 
consistent comparison of the medium effects. Our results 
yield a small but relevant variation of the one-neutron re- 
moval cross sections using the free, Pauli corrected, and 
phenomcnological NN cross sections, namely 87.62 mb, 
81.45 mb, and 81.15 mb, respectively. As observed in 
the case of the 12 C( 24 0, 23 0)X at 920 MeV/u, the use of 
Brueckner corrected NN cross sections yields 120.03 mb, 
about 37% relative to the calculations using the free NN 
cross sections. For the same reason as with the previ- 
ously considered reaction, this discrepancy is meaningless 
as one extrapolates the Brueckner results beyond their 
regime of validity. 

Both reactions considered above are very illustrative 
as they show a great consistency between the calcula- 
tions performed by different authors, with somewhat dif- 
ferent methods. They also show the expected relevance 
of medium corrections at intermediate and low energy 
collisions. 



G. Relevance of medium effects 

As we mentioned above, medium effects have been rou- 
tinely neglected in the experimental analysis of knockout 
reactions. But their relevance has been known for a long 
time in the analysis of elastic and inelastic scattering, as 
well as of total reaction cross sections [13, EH . The ef- 
fects are larger at lower bombarding energies, where Pauli 
blocking strongly reduces the nucleon-nucleon cross sec- 
tions in the medium. A systematic study of these effects 
has been presented in Ref. (l8l |. 

To corroborate these statements, in figure [TTJ we show 
the data on p + 12 C reaction cross sections taken from 
the Ref. [4l[ in the energy region of our interest, 20-100 
MeV/nucleon. The cross sections were calculated from 
the relation 



<Tjt = 2tt / db b 



l-\S(b)\' 



(14) 



where S(b) has been calculated using Eqs. (|4I5I6[) and the 
carbon matter density from a Hartree-Fock-Bogoliubov 
calculation (43|- Several distinct calculations are shown. 
The solid curve uses Eq. © with the free nucleon- 
nucleon cross sections and the carbon matter density 
from a HFB calculation , whereas the triangle-dotted 
curve (the triangles are not data, but used for better vis- 
ibility) uses a different HFB densi ty J43j , consistent with 
the calculations presented in Ref. [45| . As expected, that 
the agreement between the two calculations is very good. 

The other curves in figure [TT] show the same calcu- 
lation procedure, but including medium corrections for 
the nucleon-nucleon cross section. The results are evi- 
dently very different than the previous ones. The dot- 
ted (dashed-dotted) [dashed] curves use phenomcnolog- 
ical (Brueckner) [Pauli geometrical] recipes for medium 
effects on the cross sections. Based on the large error 
bars and spread of the experimental data, it is hard to 
judge what model adopted for medium corrections yields 
the best agreement with the data. It is clear that the 
inclusion of medium effects change the results to yield a 
closer reproduction of the data. 

These findings arc in agreement with our present un- 
derstanding of medium modifications of the reaction 
cross sections and of several other reaction channels in- 
volving heavy ion scattering at intermediate energies 
(~ 50 MeV/nucleon). 



IV. CONCLUSIONS 

Often neglected effects, such as medium modifications 
of the nucleon-nucleon cross sections and Coulomb dis- 
tortion, modify appreciably the nucleon knockout cross 
sections. As we have shown, these effects do not lead 
to an appreciable modification of the shapes of momen- 
tum distributions. This is explained by the fact that the 
momentum distributions are largely the Fourier trans- 
forms of the contributing parts of the single-particle 
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wavcfunctions, overwhelmingly their asymptotic regions, 
which are the Whittakcr functions for protons or the 
Hankcl functions for neutrons, sensitive only to the or- 
bital momentum and the nucleon binding energies. We 
have shown these features explicitly by comparing our 
results with a large number of available experimental 
data. As expected on physics grounds, these correc- 
tions are larger for experiments at lower energies, around 
50 MeV/nucleon, and for heavy targets. 

As more experiments make use of heavier targets, it 
is worthwhile to illustrate the relevance of Coulomb cor- 
rections. Medium effects in knockout reactions have also 
been frequently ignored in the past. We show that they 
also have to be included in order to obtain a better ac- 



curacy of the extracted spectroscopic factors. Although 
these conclusions might not come as a big surprise, they 
have not been properly included in many previous exper- 
imental analyses. 
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